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Abstract 

The photoproduction of Drell-Yan lepton pairs at HERA is studied. We show 
that in the backward rapidity region the cross section is strongly sensitive to the 
small- x behaviour of the quark distributions in the proton, and that with sufficient 
luminosity it should be possible to distinguish singular xq ~ x~ 1 / 2 behaviour from 
standard Regge xq ~ x° behaviour. The effect on the event rate of experimental 
acceptance cuts is studied with the help of a Monte Carlo event generator. 
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1 Introduction 



The Drell-Yan process in hadron-hadron collisions has for many years yielded important 
information on the quark distributions of hadrons [|l| . In proton- nucleon collisions, for 
example, the sea quark distribution at medium and large x - which is overwhelmed in 
deep inelastic scattering by the valence quarks - can be directly measured. It was pointed 
out several years ago that the Drell-Yan process at HERA could also be of interest. In 
particular, since the photoproduction of lepton pairs arises mainly from qq annihilation, 
the quark distributions in the photon can be studied. For example, consider the cross 
section differential in both lepton-pair mass M and rapidity y, with the proton direction 
defining y > 0. In the centre-of mass frame of the ^p collision in leading order, the 
momentum fractions of the annihilating quarks are given simply by x 7 = i/rexp(— y), 
x p — V^" ex P(2/)' where r = M 2 / ' s lp . At large negative rapidities, therefore, the quark 
distribution in the photon is probed at large x, while that in the proton is probed at small 
x. Now at HERA, the fact that the photon-proton system has a large positive boost 
in the lab frame means that yi^ < for the lepton pairs already implies a very large 
rapidity in the photon-proton centre-of- mass frame. In fact, as we shall demonstrate 
below, lepton pairs with y^ ~ —1 probe x 1 ~ 1 and x p ~ 10~ 4 — 10~ 3 . In addition, 
such lepton pairs should be relatively easy to identify, being produced in the 'backward' 
part of the detector which is relatively free from hadronic activity. 

The importance of such a measurement is that it allows an independent probe of the 
small x structure of the proton. Recent deep inelastic structure function measurements 
from the HI |J and ZEUS collaborations at HERA indicate that F2(x,Q 2 ) rises 
steeply at small x and modest Q 2 , in line with predictions based on the "perturbative 
pomeron" of QCD ||. It is vitally important to obtain additional evidence of this effect, 
and the Drell-Yan cross section at negative yi^ provides just such a measurement. We 
note that the jp Drell-Yan cross section has already been advocated as a probe of the 
small x structure of the photon 0, |6|, [7j , but this corresponds to large, positive y\^ values. 
In practice, the hadronic activity in this region presumably makes a clean measurement 
rather difficult. 

In this paper we study the jp Drell-Yan cross section for yi a b < 0. To test the dis- 
criminating power in the i p < 1 region, we use the recent MRS-Dq and MRS-D'_ parton 
distributions [|J. The former have "conventional" Regge small x behaviour, xq,xg ~ a; 
as x — > 0, while the latter have more singular xq, xg ~ x~ x l 2 behaviour as expected 
from perturbative pomeron arguments || . Note that the recent F2 data from HI || and 
ZEUS H at HERA strongly favour the D'_ distributions. A crucial part of our argument 
is that for yi^ < the parton distributions of the photon are being probed in a x region 
where they are already constrained by K7 data. The difficulty here stems from the fact 
that in practice the K7 data in the relevant x, Q 2 regions are not very precise, and this 
introduces some uncertainty into the Drell-Yan predictions which in turn weakens the 
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sensitivity to the proton structure. To try to quantify this, we use a variety of recent 
photon parton-distribution parametrizations. 

The paper is organised as follows. In the next section we describe the theoretical 
formalism for computing the Drell-Yan cross section at next-to-leading order in photon- 
proton collisions. Numerical calculations then quantify the overall event rates and the 
differences obtained from using different sets of parton distributions. In Section 3 we 
perform a Monte Carlo event simulation to study in detail the x values which are being 
probed, and also to investigate the energy and angular distributions of the leptons in the 
lab. Section 4 contains our conclusions. 



2 The Drell-Yan cross section at next-to-leading or- 
der 

2.1 Expression for the inclusive cross section 

In this study we use both next-to-leading order matrix elements and parton distributions, 
in the MS scheme. We calculate the inclusive cross sections of dilepton production at 
HERA: 
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where r = M 2 / s and the higher-order terms f qq - and f gq are given in Eqn.(5) in ||. 
In this expression, the first term contains the dominating LO qq initial state and the 
associated real and virtual corrections. The second term contains the g p + g 7 and q p + g 1 
initial states and finally the third term is the Compton term given by q p + 7 — > 7 * +q. 
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Throughout this calculation we set the factorization and renormalization scales at the 
invariant mass of the dilepton pair; that is, fi 2 = Q 2 = M 2 . 

For the parton distributions in the photon which appear in the above expression, we 
take the GRV(NLO) next-to-leading order set as standard and probe the sensitivity of 
the cross sections to the MRS-Dq and MRS-D'_ parton distributions. We also study the 
sensitivity of the cross sections to different photon parametrizations, to attempt to quan- 
tify the uncertainty in extracting information on the small-x proton distributions from 
this source. To obtain the relevant cross sections for the ep initial state we convolute the 
cross-section expression with the simplest form of the equivalent photon approximation 



TO], in which we have set the energy scale to be s = x p x e s. 



2.2 Event rates and sensitivity to photon and proton parton 
distributions 

The cross-section calculations are performed for the energy parameters of HERA. For 
the 7p initial state we take \/s = 200 GeV, and for the ep initial state we assume 
E e = 30 GeV and E p = 820 GeV. In Figure 1 we present the results of the calculation 
of the cross section dajdM as a function of M, where M is the invariant mass of the 
dilepton pair. The 7p cross section is presented in Figure 1(a) for both the MRS-Dq and 
MRS-D'_ parton distributions. The MRS-D'_ cross section is larger, as expected, but 
the difference between the two decreases substantially as M increases, since at small M 
the sensitivity to small-x in the proton is large and the differences between the two sets 
of proton distribution functions is accentuated. At M = 2 GeV/c 2 the cross section is 
approximately 76 pb/GeV for the MRS-D case and 155 pb/GeV for the MRS-D'_ case, 
and at M = 4 GeV/c 2 the cross sections have fallen to 10 pb/GeV and 14 pb/GeV, 
respectively. The next-to-leading contribution to this cross section is approximately 1.5 
to 2 times that of the leading order one. In Figure 1(b) the corresponding cross section 
for the ep initial state is shown. Evidently these cross sections are approximately one 
order of magnitude smaller than those for 7p. The reason is that at these energies 
C-y P ~ 0.1£ ep . The sensitivity of the ep cross sections to the two sets of parametrizations 
is somewhat smaller than in the 7p case. 

To emphasize the sensitivity of the Drell-Yan cross sections to the small-x behaviour 
of the proton distribution functions we next calculate the rapidity distribution of the 
dilepton pairs in the jp case. The kinematics are such that x p = y/re u ' yp and x 7 = 
■s/re~ y ' yp , which implies that x p — > in the backward direction, which is the region where 
differences between the MRS-D and MRS-D'_ should be observed. In Figure 2 we present 
the plots for the cross section da/dMdy as a function of y for M = 4, 6 GeV/c 2 . The 
M = 4 GeV/c 2 case is shown in Figure 2(a) and we see that at large negative y (= y^b), 
the MRS-D'_ cross section is larger than the MRS-D one by a factor larger than 2. At 
this invariant mass the proton is probed down to x p ~ 4 x 10~ 4 (see below). At smaller 
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M the difference is considerably larger; for instance at M = 2 GeV/c 2 there is a factor of 
4 between the two distributions at large negative rapidities. However, at such low masses 
data from hadron-hadron collisions suggest that the background to the dilepton cross 
section arising from J/ip production may be substantial 0. In what follows, therefore, 
we impose a lower limit of M > 4 GeV/c 2 on the dilepton mass. In Figure 2(b) the same 
distribution is plotted for M = 6 GeV/c 2 and as expected, the effect washes out, due to 
the fact that a larger invariant mass forces larger x p values; at large negative rapidities 
the cross section ratio, MRS-D'_/MRS-Dq reduces to approximately 1.5. In the positive 
rapidity region the ratio is very nearly 1 in both Figures 2(a),(b). 

The usefulness of this process to probe the small-x behaviour of the proton structure 
function depends on the insensitivity of the cross sections to the large- x behaviour of the 
photon structure function. To investigate this, we repeat the calculations using several 
other sets of photon distributions, keeping the proton distributions (MRS-D'_) fixed. 
Figure 3 shows the rapidity distribution, da/dMdy for M = 4 GeV/c 2 as a function of y 
for various choices of the photon distributions: (i) the reference GRV(NLO) set [[| (solid 
line), (ii) the GS(NLO) set (dashed line) and the LACl(LO) set @ (dotted line). In 
the interesting region of phase space (large negative rapidities) we see that the different 
distributions give roughly the same shape in rapidity but a spread in normalisation of 
order 10%. This can be understood in terms of the spread in predictions of these sets 
for the photon structure function F%. Figure 4 shows data on F% together with the 
predictions of the same three sets as in Figure 3. We have chosen data in a similar 
Q 2 range as the M 2 values relevant to the Drell-Yan cross section at HERA. We see 
immediately that the spread in the predictions simply reflects the uncertainty in the 
structure function measurements. More precise measurements of F% would of course pin 
down the Drell-Yan cross section more accurately. In the meantime we should regard 
the 0(10%) spread in normalisation as a systematic uncertainty in the extraction of the 
proton distributions from the Drell-Yan data. 



3 Event simulation 

In the previous section we have seen how the jp Drell-Yan cross section for y lab < 
and M > 4 GeV/c 2 offers a good possibility of obtaining information on the quark 
distribution in the proton at small x. An important question is whether the bulk of 
these events actually fall within the acceptance of the HERA detectors. To study this, 
we perform a Monte Carlo simulation using the PYTHIA5.6 [[15] MC generator for ^p — > 
l + l~ + X at y/s~ip = 200 GeV. This simulation also allows us to investigate which x p , x 7 
values give the dominant contribution to the cross section. 

Figure 5 show the distributions in (a) x 7 and (b) x p for M > 4 GeV/c 2 . The two 
histograms correspond to < (solid lines) and yi a b > (dashed lines). The former 
is dominated by x p values in the 10 -4 — 10 -3 region, as expected. Selecting only those 
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events with M > 4 GeV/c 2 , and denning $1 and E t to be the lepton polar angle and 
energy in the lab frame respectively, Figure 6 shows the scatter plots for (a) x^ vs. $1, 
(b) x p vs. (c) Ei vs. $1, and (d) M vs. We see immediately that most of the 
events have leptons which are produced with reasonable energy and at sizeable angles to 
the beam direction d = 0°, 180°. Given that we expect very little additional hadronic 
activity in these regions of phase space, the acceptance for both muon and electron pairs 
should be quite high. 

4 Conclusions 

The structure of the proton at small x is an important new area of interest, and the first 
F 2 data from HERA provide tentative evidence for the "perturbative pomeron". The 
Drell-Yan cross section in the backward rapidity region should in the very near future 
give a confirmation of this small- x behaviour from a completely different process. The 
decrease in hadronic activity as y^h becomes more negative and the effect of the boost on 
the lepton lab angle should make this a particularly clean measurement. Our study has 
concentrated on 7^ collisions at v /s^p = 200 GeV, tacitly assuming that the final-state 
electron can be tagged with reasonable efficiency at HERA. However, as is clear from 
Figure 1(b), the effects we describe should also be observable in the ep cross section. 

Observation of the dramatic differences induced by "standard" and "singular" small- 
x behaviour in the proton in the yi a b < Drell-Yan cross section relies on knowledge 
of the large-x photon structure. We have shown that the uncertainty in the photon at 
large x is of the order of 10%, which for M = 4 — 6 GeV is small in comparison to the 
Dq/D'_ effect. But we should stress again that once these Ff become well-measured at 
HERA, the Drell-Yan cross section will provide a complementary measurement of the 
large-x behaviour of the photon. 
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Figure Captions 



[1] Predictions for the Drell-Yan cross section da/dM in (a) 7p collisions at y/s lp = 200 
GeV, and (b) ep collisions at v /s^ = 314 GeV. The curves correspond to next-to- 
leading order calculations using the GRV(NLO) distributions for the photon || 
with MRS-Dq (dashed line) and MRS-D'„ (solid line) distributions for the proton 

!■ 

[2] Predictions for the Drell-Yan rapidity distribution da/dMdy for (a) M — 4 GeV/c 2 
and (b) M = 6 GeV/c 2 in jp collisions at yfs^ = 200 GeV. The curves correspond 
to next-to-leading order calculations using the GRV(NLO) distributions for the 
photon H with MRS-D' (dashed line) and MRS-D'_ (solid line) distributions for 
the proton ||. 

[3] Predictions for the Drell-Yan rapidity distribution da/dMdy for M = 4 GeV/c 2 in 
7P collisions at ^Js lp = 200 GeV. The curves correspond to next-to-leading order 
calculations using (a) GRV(NLO) @ (solid line), (b) GS(NLO) (dashed line) 
and (c) LACl(LO) @ (dotted line) distributions for the photon, with MRS-D'_ 
distributions for the proton M. 



[4] The photon structure function F%{x, Q ). The data are from the PLUTO ]T^] and 



TASSO 0] collaborations. The curves are the theoretical predictions using the 
same sets as in Figure 3: (a) GRV(NLO) § (solid line), (b) GS(NLO) (dashed 
line) and (c) LACl(LO) Q (dotted line). 

[5] Distributions in (a) x 7 and (b) x p for M > 4 GeV/c 2 , using a PYTHIA5.6 |J MC 
simulation of ^p — > l + l~ + X at y/s 1P = 200 GeV. The two histograms correspond 
to 2/iab < (solid lines) and yi^ > (dashed lines). 

[6] Scatter plots for (a) x 7 vs. (b) x p vs. $i, (c) Ei vs. and (d) M vs. 
for M > 4 GeV/c 2 , using a PYTHIA5.6 pf MC simulation of jp -> Z+r + X at 
,/sZ:= 200 GeV. 
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